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Abstract. We briefly summarize our present knowledge of the theory of low-mass 
stars and substellar objects and their contribution to the Galactic population. 



1 Introduction 

The search for substellar objects (SSO) has bloomed over the past few years 
with the unambiguous identification of several free floating brown dwarfs 
(BDs) and with the discovery of numerous planets orbiting stars outside the 
solar system. Substellar objects therefore exist and ongoing and future obser- 
vational projects are likely to reveal dozens more of such objects. Two groups, 
the Lyon group and the Tucson group, have aspired to derive a complete the- 
ory of the evolution and the spectral signature of low-mass, dense objects, 
from Sun- like to Saturn- like masses, covering 3 orders of magnitude in mass 
and 9 in luminosity, bridging the gap between stars and gaseous planets (see 
e.g. Burrows et ai, 1997 and Baraffe et ai, 1998 and references therein). 
These two groups have incorporated the best possible physics aimed at de- 
scribing the mechanical and thermal properties of these objects - equation 
of state, synthetic spectra, non-grey atmosphere models. Both groups have 
"met in the middle" by successfully identifying the main spectral properties 
of the benchmark BD G1229B, providing a determination of its mass and age 
(Marley et ai, 1996; AUard et ai, 1996). We refer to the afore-mentioned 
papers and the references therein for a complete description of the physics 
entering these models. Only a brief outline of these characteristics is given 
below. 



2 The physics of sub-stellar objects 
2.1 Interior 

Central conditions for SSOs are typically Tc^lO^ K and pc ^ 10^-10^ g 
cm~^, characterizing a strongly coupled electron- ion plasma. Moreover, in 
the envelope, the electron binding energy is of the order of the Fermi energy 
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Ze^ /ttQ ^ Ep so that presswre-dissociation and ionization take place along 
the interior. Very recently laser-driven shock wave experiments at Livermore 
have achieved pressures on liquid D2 of up to 5 Mbar at high temperature 
(Collins et al., 1998), exploring for the first time the regime of pressure- 
dissociation and ionization, therefore probing the equation of state (EOS) 
under conditions characteristic of SSO interiors. As shown in Figure 3 of 
Collins et al., the experimental hugoniot revealed the excellent behaviour of 
the EOS developed by Saumon and Chabrier (Saumon et al., 1995 and refer- 
ences therein). In particular the pronounced compressibility observed in the 
dissociation domain (p/po = 5.88) agrees remarkably well with the theoreti- 
cally predicted value. These experiments open a new window in physics and 
astrophysics by constraining the physics of the interior of SSOs in laboratory 
experiments. 

2.2 Atmosphere 

In the atmosphere, the hydrostatic equilibrium condition dr = RdP/ g, where 
g = Gm/R^ ~ lO'^-lO^ cm s~^ is the characteristic surface gravity of SSOs, 
yields Pph ^ g/n «0.1-10 bar and pph ~10~^-10~^ g cm~^ at the photo- 
sphere. Collisional effects are important under such conditions and introduce 
sources of absorption like e.g. the collision-induced absorption (CIA) between 
H2-H2 or H2-He below ~ 5000 K (see e.g. Borysow et al., 1985). In the ef- 
fective temperature range characteristic of low-mass stars (LMS) and SSOs 
(Teff;^5000 K), numerous molecules form, in particular metal oxydes and 
hydrides (TiO,VO,FeH,CaH), the major absorbers in the optical, and CIA 
H2, H2O, CO which dominate in the infrared (see Allard et al., 1997 for a 
review). The situation becomes even more complicated for SSOs, due to the 
changes in molecular chemistry which occur across the atmospheric temper- 
ature range from 2000 K (the coolest stars) to 170 K (jovian conditions). 
At 2000 K, most of the carbon is locked into carbon monoxide CO, while 
the oxygen is found in water vapor H2O, dominantly, and in titanium TiO 
and vanadium VO monoxides. Below ^ 1800 K, the dominant equilibrium 
form of carbon is no longer CO but CH4 (Fegley & Lodders, 1996). As con- 
firmed by the observation of G1229B (Oppenheimer et al., 1995), methane 
features begin to appc;ar in the infrared while titanium dioxide and silicate 
clouds form at the expense of TiO, modifying profoundly the opacity of the 
atmosphere. For jovian-like atmospheres, the dominant equilibrium form of 
nitrogen is NH3 {Tea < 600 K) and below Tdi ~ 200 K water condenses to 
clouds at or above the photosphere. As shown by Tsuji et aL(1996) and Jones 
& Tsuji (1997), there is evidence for condensation of metals and silicates into 
grains (e.g. TiO into CaTiOs, Mg, Si into MgSiOa) for Teff;^3000 K, i.e. at 
the bottom of the main sequence. 

A proper inclusion of this complex atmospheric chemistry, and of a proper 
calculation of unknown cross-sections, represent the main challenge for a cor- 
rect description of the spectral signature and the evolution of substellar ob- 
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jects. This grain formation process has been included in the calculations of 
the Tucson group only in term of condensation. The species are precipitated 
according to their condensation curves. If a species has condensed, it is left 
at its saturated vapor pressure (Burrows et al, 1997). There is no inclusion 
of the radiative transfer effect of the grains. These models can thus be con- 
sidered as grainfree opacity models. Recently, the Lyon group has extended 
similar calculations by implicitely including all condensed species into the 
radiative transfer equations (see AUard & Hauschildt, these proceedings). 
The atmosphere profiles were matched with the interior profiles at an optical 
depth deep enough to lie on the internal adiabat (see Chabrier & Baraffe, 
1997). 

3 Color-magnitude diagram 

The present evolutionary calculations have been conducted with three differ- 
ent atmosphere models: (i) for hot objects, i.e. massive or young enough, to 
preclude the formation of grains (Teg ^ 2800 K), we have used the most recent 
grainless NGEN atmosphere models (Hauschildt et al., 1999); (ii) for objects 
in the range 3000-1000 K, i.e. from the bottom of the MS down to Gliese229B- 
like objects, we have used complete atmosphere models wich include the grain 
opacity sources in the transfer equations, the so-called DUSTY models (Al- 
lard & Hauschildt, these proceedings); (iii) for objects below 3000 K down to 
jovian temperatures, we have also considered cases where the condensates set- 
tle rapidly below the photosphere and - although modifying the atmosphere 
EOS - do not participate to the opacity, the so-called COND models. This is 
similar to the Burrows et aL(1997) calculations and is motivated by the rela- 
tive absence of grain features in the atmosphere of objects below Teg ~ 1000 
K, i.e. Glicsc229B-likc and cooler objcx'ts. We fomid the cfFcx-t of grain opac- 
ity (DUSTY) to affect only moderately the H-burning minimum mass. Mod- 
els with grainless atmosphere yield m = 0.072 M©, L = 5 x 10~^Lq and 
Teff = 1 700 K at the H-burning limit, whereas models with grain opacity 
give m O.O7M0, L « 4 x 10"^ L© and T^s « 1600 K, for solar composi- 
tion. 

Figure 1 displays a K vs ,J-K color-magnitude diagram (CMD). In terms 
of colors there is a competing effect between grain and molecular opacity 
sources for objects at the bottom and just below the MS. As already identified 
for e.g. G1229B (Allard et al., 1996; Marley et al., 1996) CH4 and CIA of 
H2 lead to bluer infrared colors for SSOs, as illustrated by the NEG and 
COND models, whereas grain opacity results in a severe redening. Indeed, 
the present DUSTY models reproduce the spectra of the DENIS objects 
(Tinney et al., 1998) and of the long puzzling object GD165B, which lies 
at the very edge of the H-burning limit (Kirkpatrick et al, 1999). As T^g 
decreases, grains settle below the photosphere and the DUSTY track will 
merge with the COND tracks. Observations of objects in the region between 
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Fig. 1. Ma' vs (J-K) diagram for 3 isochrones, namely 5x10 yrs, 5x10 yrs and 
1.2x10^ yrs (the age of the Pleiades). The small circles and squares correspond 
to the observation of main sequence stars (Leggett, 1996) and of Pleiades objects 
(Bouvier et at, 1998), respectively. Some identified BDs are also indicated 



these two extreme cases will be the next confirmation of the present theory. In 
spite of this strong absorption in the IR, however, BDs around 1500 K radiate 
nearly 90% (99% with dust) of their energy at wavelengths longward of 1 iim 
and infrared colors are still preferred to optical colors (at least for solar metal 
abundance), with M, L, K as the favored bands, and Mm ~Ml ~ 11, ~ 
12 at the H-burning limit, at 5 Gyr. Direct observation of the characteristic 
fluxes of SSO's are now within reach with several observational projects, like 
SOFIA, SIRTF, ISAAC (Burrows et al, 1997; AUard & Hauschildt, these 
proceedings). 

Several BD surveys are presently conducted in young clusters and it is im- 
portant to develop accurate (non-grey) models for pre-MS stars and young 
BD's. The more massive of such objects will be hot enough so that grain 
formation does not occur, but for t<bx 10^ yrs, proper evolutionary calcu- 
lations must consider two effects : (i) first of all the influence of the initial 
conditions must be carefully examined because they will significantly affect 
the mass-magnitude relationship; (ii) for young objects, the gravity is small 
(log (7 ^3) and sphericity effects in the resolution of the transfer equations 
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might come into play. The derivation of such complete models is under work, 
but Figure 2 displays the present pre-MS models (Baraffe et al., 1998) based 
on non-grey plane-parallel atmosphere models in a theoretical HR diagram. 




Fig. 2. Pre-main sequence isochrones (models are available from ftp site - see at 
the end of the paper) 



4 Stellar and brown dwarf mass functions 
4.1 Stellar mass function 

The determination of the contribution of LMS and SSOs to the Galactic mass 
budget requires the correct determination (slope and normalization) of the 
mass- function (MF) down to the hydrogen burning limit, in order to have 
a solid foundation for extrapolation into the substellar domain. This issue, 
however, is presently not completely settled for the Galactic disk, and signif- 
icant differences exist between the MF inferred from parallax-determined 
luminosity functions (LF) (Kroupa, 1995) and from the HST photomet- 
ric LF (Gould, Flynn & Bahcall, 1997) (see Mera, Chabrier & Schaeffer, 
1998, Fig. 1). It is relatively safe to say, however, that the MF keeps ris- 
ing down to the H-burning limit, although with a slope shallower than a 



6 G. Chabrier et al. 



Salpeter MF. The LMS (m < 0.6 M©) mass densities inferred from the inte- 
gration of the two afore- mentioned MFs (nearby and HST) yield, respectively, 
PLMS ~ (1.9 ±0.1) X 10-2 and plms ~ (2.6 ±0.2) x lO'^ Mq^c'^ (Mera et 
al., 1998). Adding up the contribution from more massive stars and stellar 
remnants yields « (4.0-4.6) ± 0.3 x 10"^ M0pc~^ in the Galactic disk, 
i.e. a surface density 17^, « 31 ± 2 MqPc"^. This corresponds to a stellar 
ntimfeer-density ~ 0.06 (HST) and ^ 0.35 pc""^ (nearby). 

For the spheroid the question is more settled with a MF with a<^l (where 
dN /dm (X m~") and a stellar density < 4.0x 10~^ Mqpc~^, less than 1% of 
the required dynamical density (Graff & Freese, 1996; Chabrier & Mera, 1997; 
Gould et al., 1998), i.e. a number density < 10""^ pc""^. Extrapolating into 
the BD domain yields for the spheroid ubd < 10~^ pc~'^, pbd < /5*/10, i.e. a 
microlensing optical depth r ~ 10~^, about 1% of the value measured toward 
the LMC. The puzzle remains unsolved for the dark halo MF, although both 
the HDF observations and the narrow-range of the observed time-distribution 
of the microlensing events towards the LMC strongly suggest an IMF different 
from a Salpeter one below 1 Mq (see e.g. Chabrier, 1999). 

4.2 Brown dwarf mass function 

A proper census of the number of brown dwarfs (m < 0.07-0.09 Mq, depend- 
ing on the metallicity) has significant implications for our understanding of 
how stars and planets form. The determination of the BD MF is a compli- 
cated task. By definition BDs never reach thermal equilibrium and most of the 
BD's formed at the early stages of the Galaxy will have fainted to very low- 
luminosities {L oc 1/t, see e.g. Burrows & Liebert, 1993). Thus observations 
will likely be biased towards young and massive BD's. The age uncertainty 
is circumvented when looking for the BD MF in stellar clusters since objects 
in that case are likely to be coeval. The Pleiades cluster has been extensively 
surveyed and several BDs have been identified down to ~ 25 Mj (Martin et 
al, 1998; Bouvier et al., 1998). A single power-law function from ~ 0.4 to 
0.04 Mq seem to adequately reproduce the observations with some remaining 
uncertainty in the exponent a ~ 0.6 — 1.0 (Bouvier et at, 1998; Martin et 
al., 1998). By examinating the various observations of Doppler radial velocity 
surveys, Basri and Marcy (1997) reach the conclusion that all data are con- 
sistent with ;^ a 1 (see Basri, these proceedings for an updated analysis of 
such data). On the other hand Mayor et aZ.(1997) find for the mass function 
of the companions of the G and K dwarf sample observed with CORAVEL a 
power-law function with a ~ 0.4 from 0.4 down to 0.005 Mq. The determina- 
tion of the MF in young clusters is potentially very interesting. Indeed young 
clusters did not have time to experience evaporation in the outer regions or 
mass segregation in the central regions and the present-day MF should re- 
flect relatively closely the initial MF. Such a MF determination, however, is 
hampered by two observational and theoretical problems. From the obser- 
vational point of view, young objects are immersed into dust and a proper 
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determination of the LF requires a correct determination of the differential 
reddening in the cluster. From the theoretical point of view, none of the MFs 
determined up to now in young clusters include a careful examination of the 
effects mentioned in §3 (initial conditions, sphericity, non-grey effects) and 
thus are of dubious validity. It is thus certainly premature to claim a robust 
determination of the MF in young clusters. 

An independent, powerful information on the stellar and substellar MF 
comes from microlensing observations. Indeed, the time distribution of the 
events provides a (although model-dependent) determination of the mass 
distribution and thus of the minimum mass of the dark objects: dNev/dte = 
E X e{te) X dr/dte oc P{m) / y/m, where E is the observed exposure, i.e. the 
number of star X years, e is the experimental efficiency, F is the event rate and 
P{m) is the mass probability distribution. The analysis of the published 40 
MACHO + 9 OGLE events towards the bulge is consistent with a rising mass 
function at the bottom of the MS with a minimum mass rriinf ^ 0.05 M©, 
whereas a decreasing MF below 0.2 Mq is excluded at the 95% confidence 
level (Han & Gould, 1995; Mera et al., 1998). Although the time distribution 
might be affected by various biases (e.g. blending) and robust conclusions 
must await for larger statistics, the present results suggest that in order to 
explain both star counts and the microlensing experiments, a substantial 
amount of SSOs must be present in the galactic disk. Indeed, extrapolation 
of the stellar MF (§4.1) into the BD domain down to 0.05 Mq yields for the 
Galactic disk pbd ^ 4.0 x lO"''' Mgpc"^, i.e. Sbd ~ SMgpc'^, i.e. a BD 
number-density comparable to the stellar one, ubd O.lpc"^ ~ n*. 

For the spheroid, extrapolation of the afore-mentioned stellar MF yields 
PBD < 10~^ Mqpc"^, nsD < 10~^pc~^, whereas for the dark halo the nor- 
malization is about 2 orders of magnitude smaller (Chabrier & Mera, 1997). 

4.3 Planet mass function 

It is obviously very premature to try to infer the mass distribution of ex- 
oplanets. This will first require a clear, both theoretical and observational, 
distinction between planets and brown dwarfs. An interesting preliminary re- 
sult, however, comes from the observed mass distribution of the secondaries 
conducted by Mayor and collaborators. As shown by Mayor et aZ.(1997), there 
is a strong discontinuity in the mass distribution at m2/sini « 5Mj, with 
a clear peak below this limit. If confirmed this would suggest that planet 
formation in a protoplanetary disk is a much more efficient mechanism than 
BD formation, which results from cloud collapse and fragmentation. 

5 Conclusion and perspectives 

As mentioned in the introduction, accurate models for low-mass stars, brown 
dwarfs and giant planets are needed to shed light on the observable properties 
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of these objects and to provide guidance to the ongoing and future surveys 
aimed at reveaUng their contribution to the Galactic population. From this 
point of view, tremendous progress has been made in the recent years with 
the derivation of consistent evolutionary model and synthetic spectra calcula- 
tions which accurately reproduce the observed sequences of globular clusters, 
Pleiades, field objects and of the benchmark BD G1229B in various photomet- 
ric passbands (see references mentioned in the Introduction). Moreover the 
afore-mentioned LMS models have been shown to accurately reproduce the 
observational mass-magnitude relationship (Henry & McCarthy, 1993; Henry 
et al., 1999) both in the optical and in the infrared (Baraffe et al, 1998) 
and to yield a consistent, coeval sequence for the quadruple system GG-Tau, 
whose component masses extend from 1.2 Mq down to ~ 0.035 Mq (White et 
al, 1999), a formidable test for the theory. On the other hand, stringent con- 
straints on the theory of dense/cool objects are now provided by laboratory 
high-pressure experiments, for the interior, and by various spectroscopic and 
photometric observations of LMS and SSO's. Any theory aimed at describing 
the mechanical and thermal properties of these objects must be confronted 
to these experimental/observational constraints in order to assess any degree 
of validity. Improvement in the theory will proceed along with the discovery 
of many more SSO's, hopefully bridging the gaps on either side of G1229B 
from the bottom of the MS to Jupiter-like objects. Observation of a transit 
EGP would allow the determination of the radius and the mass of the object, 
providing a formidable constraint on the theory. For 51-Pcg-likc objects, the 
probability to observe such a transit p = di,/D, where and D arc the stellar 
and orbital diameters, respectively, is p ~ 10%, by no means negligible. 

The increasing number of observed LMS and SSOs, together with the 
derivation of accurate models, will allow eventually a robust determination 
of the stellar and substellar mass functions, and thus of the exact density of 
these objects in the Galaxy. As discussed in §4, present MF determinations 
in various Galactic regions point to a slowly rising MF near and below the 
H-burning limit, with a BD number density comparable to the stellar one. 
Whether this behaviour is universal (although we already know it is certainly 
not the case for the dark halo), whether it is consistent with a general log- 
normal form, must await confirmation from future observations. On the other 
hand, the amazingly rapid pace of exoplanet discoveries should yield the 
determination of the planetary MF. These combined informations will allow 
the determination of the BD minimum mass and planet max;imum mass. 

As seen, the physics of SSOs involves an amazingly large domain of physics 
and astrophysics, from the fundamental N-body problem to star formation 
and Galactic evolution, and will certainly remain a lively domain of the early 
21st century astronomy. 

The present grainless models are available from: 

ftp ftp.ens-lyon.fr, username: anonymous 
ftp > cd /pub/users/CRAL/ibaraffe 
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ftp > get BCAH98_models and BCAH98.inodels_BD 
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